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Abstract 

A discrete version of the inverse scattering method proposed by Ablowitz and Ladik is 
generalized to study an integrable full-discretization (discrete time and discrete space) of 
the coupled nonlinear Schrodinger equations. The generalization enables one to solve the 
initial-value problem. Soliton solutions and conserved quantities of the full-discrete system 
are constructed. 

1 Introduction 

Among various soliton equations, a system of coupled nonlinear Schrodinger (CNLS) equations, 



■ d 9j . d 2 qj v-^ 



. dr<i d 2 Vj -l— 



fc=i 



j = l,2,...,m, (1) 



dt dx 



k=l 



is physically significant since it describes multi-mode wave propagation in nonlinear optics. The 
CNLS equations can be solved via the inverse scattering method (ISM) under the reduction rj = 
—q* and the decaying boundary conditions qj — > as x — > ±oo. It is, however, not always 
easy to trace time evolution of initial values for soliton equations solvable via the ISM. Thus, 
it is an important problem to hnd good schemes for numerical computation of soliton equations. 
One possible approach is to find discretization of soliton equations with preserving the complete 
integrability (integrable discretization for short). The discrctized soliton equations are not only 
useful for numerical computation but also interesting as models of nonlinear population dynamics, 
nonlinear electric circuits and nonlinear lattice vibrations. 

An integrable discretization of the nonlinear Schrodinger (NLS) equation ((Q) with m — 1) was 
proposed by Ablowitz and Ladik J2|, |). The author, Ujino and Wadati generalized the formulation 
of Ablowitz and Ladik and found an integrable semi-discretization (space discretization) of the 

* Supported by JSPS Research Fellowships for Young Scientists. 



1 



coupled modified KdV equations J4| and the CNLS equations 0. In this article, as a continuation 
of H , we consider an integrable discretization of both space and time (full-discretization for short) 
for the CNLS equations (Q). In accordance with the page limit, we concisely explain an outline 
of the ISM for the model. The reader can refer to (§ § for more details because we can perform 
the ISM in the same manner as in these papers, apart from the difference in time dependences of 
scattering data. 



2 Lax pair 

We introduce a set of auxiliary linear equations, 

* n+ i = L n f n , * n = V n V n , (2) 

where the tilde (~) denotes the time shift in discrete time I € Z (I — » I + 1). The compatibility 
condition of (^) is given by 

L n V n = V n +\L n . (3) 

L„ and V n , and (|^) are respectively the Lax pair and the zero-curvature condition (or Lax equation) 
in the full-discrete case. To obtain a full-discrete version of the CNLS equations, we choose the 
form of the Lax pair as 
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Here z is the spectral parameter which is time- independent. Q n and lZ n are p x p matrices. H\ 
and H% are p x p constant matrices, a and (3 are constants. j n is a scalar variable. We assume the 
following relations, 

[H 1 ,Q n 1l n ] = [H 2 ,n n Q n ] = 0, (6) 
Hi Q n - Q n H 2 = -2Fx Q n F 2 , H 2 TZ n - K n H 1 = 2F 2 TZ n F 1 , (7) 
where [ ■ , • ] denotes the commutator and F\ , F 2 are p x p constant matrices which satisfy 

(Fi) 2 = {F 2 f = I, [F^Hx] - [F 2 ,H 2 ] = O. (8) 

Here / is the identity matrix. We notice that the following relations hold: 

e~ ivHl Q n e iyH2 = cos(2y)Q n + ism(2y)F 1 Q n F 2 , (9a) 



e lyH2 K n e-' lyHl = cos{2y)1l n +i&w{2y)Ftfl n F 1 . 



(9b) 
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We write the difference interval of time by St. Putting (|j) and ^) into (|J), we obtain 



. 1 ~ tan (3 cot(2a + 0) , - . . s „ 

lT-(Q„ - Qn) ~ 1 fT {jn+lQn+l{I ~ UnQn) ~ Jn+l (i - QnF„ Qn-l} 

or or 



tan /3 
<5r 



{jn+l^lQn+lW - K„Q n ) +jn+l( I ~ QnTl n )F 1 Q n - 1 F 2 



F 1 (Q n + Q n )F 2 } = 0, 



(10a) 
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{jn+lF^Hn+lFiil - Q n Kn)+3n+l(I - K n Q n )F 2 K n - 1 F 1 



-F 2 {ji n + n n )F 1 } = 0, 

j n +l(I ~ QrJt n ) = jn{I ~ Qn^n), 
j n+ l{I ~ KnQn) = jn(I ~ KnQn), 

with the help of (g). If we set 

7T tan 8 



or equivalently, 



7T 1 

a = tan -1 St, 8 = tan -1 St, 

4 2 ' H 



(10a) and (10tj) are simplified as 
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(lOd) 



i^-CQn - Qn) + jn+iFiQ n+ iF 2 (I - K n Q n ) + j n+1 (I - Q„7^)F 1 Q n _ 1 F 2 
ot 

-F 1 (Q n + Q n )F 2 = 0, 



(11a) 



i-(^„ - F„) - j n+1 F 2 K n+1 F 1 (I - Q n Kn) - jn+i(I - T^nQn)F 2 n n -iFi 
ot 

+ F 2 (K n + K n )F 1 =0. 



(lib) 



To consider a reduction of these matrix equations to full-discrete CNLS equations, we define Fi, 
F 2 , Hi, H 2 , Q n and lZ n recursively as follows: 

F W - i F W - i _ _i w-(i) _ i n(i) - (i) 7?(i) - r (i) 
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Here J 2 m-i is the 2 m 1 x 2 m 1 identity matrix. We can easily show that (^))-((|) and the relation, 

m 



k=l 



are satisfied. Then it is straightforward to prove that substitution of F[ m \ F^ m} , Q ( n n} and TZ^ 
into (10c), (lOd), (11a) and (lib) yields the following full-discrete CNLS equations, 

1 m 

- #>) + (i-E aS'VW)^!^ + - (#> + <#>) = o, 

fc=l 

1 m 

i^(#> - r») - (l - E ^9^) jn+l(f »+ 1 + r " 2l) + (f " } + )} = °' 



k=l 
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m m 

(i-E = i- ( x - E aW) • i = i, 2, . . . , m. 



(12) 



fc=i fe=i 
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Assuming the reduction r„ = —q n , we obtain the self-focusing case of the full-discrete CNLS 
equations: 

1 m 

^jM ] - <#>) + (i + E \in ] l a )jn+i(gSi + - + g«>) = 0, 
k=l 



j n+1 (l + Ek fe) | 2 )=Jn(l+E^' 



s(fe)|2 
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fe=i 



In the one-component case (to = 1), this system coincides with the full-discrete NLS equation 
obtained by Ablowitz and Ladik H. If we take the continuum limit of time (St — > 0), we can 
equate j„ to 1. Thus, the system rtl2) reduces to the semi-discrete CNLS equations H, 
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fe=i 
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The corresponding linear problem is given by taking the continuum limit: 
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It should be noted that the above Lax pair for the semi-discrete CNLS equations is equivalent to 
the one given in || via a gauge transformation. 



3 Inverse scattering method 

As was shown in the previous section, the linear problem for the full-discrete CNLS equations ([l2] 
is expressed in terms of a 2 m x 2 m matrix as 
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where a = ? — i tan 1 5t. To simplify the analysis, we introduce a gauge transformation (I: 
discrete time), 



(3 = tan 1 St, 







" e -i(na+lf3)H[ m) 






$2n _ 




e -i(na+;/3)_ff<"* ) 
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by which the linear problem (114) is changed into the standard form, 
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Here the transformed potentials are given by 



If we set 



Q{m) = e -i(nQ+!/3)<» q(tu) ^(na+lfj)^" 1 
R( m ) = e -i("«+i/3)^ m) -^(m) e i( nQ +;/5)H;" 
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j = 1,2, . . . ,m, 



Qn and i?n for to > 2 are expressed as 
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Here 2 m 1 x 2 m 1 matrices {ei, • • • , e2 m -i} satisfy the following important relations: 



{e»>ej}+ 



-2% I, 



(17) 



where {■, ■}+ denotes the anticommutator and 1 is the 2 m_1 x 2 m_1 unit matrix. 

It should be remarked that the transformation (|l6|) connects the full-discrete CNLS equations 
( |l2] ) with full-discrete coupled modified KdV equations || ||, 
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where 



sin [3 



We can investigate the direct and the inverse scattering problems associated with ( |l5| ) in the 
same way as in the semi-discrete case [|], || by assuming m > 2, the reduction and 
the boundary conditions, 



7 0') r U) n i 

in i ' n u > J n 



St 
sin/3 



as n — > ±oo. 



(19) 



In fact, j„ - * const, (n — > ±oo) is sufficient to perform the ISM. We have chosen the constant as 
@ just for convenience. 

We introduce Jost functions <j> n and </>„ as solutions of (|l5|) which satisfy the boundary condi- 
tions: 
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O 
-I 



as n 



A(z)z n 




B{z)z~ n 
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£(z)z n 
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Here {A(z), -B(z), B(z)} are n-independent 2 m ~ 1 x 2 m_1 matrices which are called scattering 

data. Suppose that 1/ det A(z) has 2N isolated simple poles {zi, z 2 , ■ ■ ■ , z 2 n} in z| > 1 an d is 
regular on the unit circle C. It will be explained why we choose the number of poles to be even. 
Let us define F(n) in terms of the scattering data by 

2N 

F(n) = —f B(z)A(z)- 1 z- n - 1 dz + J2C j zr n -\ 

where Cj is the residue matrix of B(z)A(z) at z — Zj. Similarly, we suppose that 1/ det A(z) 
has 27V isolated simple poles {z\, z 2l ■ ■ ■ , z 2 n} m \ z \ < 1 an d i s regular on the unit circle C. We 
define F(n) by 

2JV 



F(n) = ±-<f S(z)A(z)- 1 z- 1 dz - £ C fc *T 1 

71-1 ^ C k=l 



where Ck is the residue matrix of B(z)A(z) 1 at z — Z}~. It can be shown after some computation 
that F(n) and F(n) are connected with the physical variables Qn™^ and ' through the following 
Gel'fand-Levitan-Marchenko equations: 

oo oo 

K 1 (n,q)=F(n + q)- ^ n ") F ( n " + n')F(n' + q), (20) 



n'=n+l n"— n+l 



K 2 (n, q) 
for q > n. Here 



-,F(n + g) — K 2 (n,n")F(n" + n')F(n + q), 

n' — n+1 n ;/ — n+1 

-Ki(n,n+ 1) = Q T ( r } , -R 2 (n,n + 1) = ^ m) . 



(21) 



It should be noted that the scattering problem ([l5]) gives the symmetry properties of the 
scattering data. Iterating (lq), we can prove that A(z), A(z) are polynomials in z of even degree 



G 



and B(z), B(z) are polynomials in z of odd degree. This fact gives the symmetry properties of F 
and F: 

2Fn(n), n : odd, 
O, n : even, 



F(n) 
F(n) 



2Fn(n), n : odd, 
O, n : even. 



Further, due to the internal symmetries of Q„ and given by (17), Fr(ji) and Fr{ti) for odd 
n are expressed as 

2m-l 2m-l 

F«(n) = / (0) (n)l+ £ / W ("H, F R (n) = (n)l - £ / (fc) (n)e fe . 

fc=i fc=i 

Here the functions f^ '(n) and are real. From the asymptotic form of the Lax matrix V„ 

as rt — » ±00, the time dependences of the scattering data are given by 

A(z,l)=A(z,0), (22a) 



/ 1 _ J)f 7 2 1 

B(z, l)A(z, I)' 1 = B(z 7 0)A(z, 0)- 1 ( 1 _ ^ j , (22b) 



/ 1 - tffc z? , 

C s (l) = C s (Q) - g-A . (22c) 



1 — tft/ 

Thus f(°\n) and f^ k \n) satisfy the following dispersion relation: 

i{/«(2n + 1) - /«(2n + 1)} = /«(2n + 3) - /«(2n - 1), j = 0, 1, . 



,2m - 1. 



We can solve the initial- value problem of the full-discrete CNLS equations ( |l3| ) in the following 
steps. 

(a) For given potentials at I = 0, Qn"^(0) and R^fi), we solve the direct problem of scattering, 
(|l5|), and obtain the scattering data {B(z)A(z)~ 1 , Zj, Cj}. 



(b) The time dependences of the scattering data are given by (|2 

(c) We substitute the time-dependent scattering data into the Gel'fand-Levitan-Marchenko equa- 
tions (20) and (pl|). Solving the equations, we reconstruct the time-dependent potentials, 

Qi m \l) = - Kl (n, n + 1: 1), R^\l) = -R 2 (n, n+1; I). 

This step corresponds to solving the inverse problem of scattering. 

In order to obtain soliton solutions, we assume that each soliton seen in J2j \Qn {l)\ 2 has a time- 
independent shape. By calculating an asymptotic behavior of the tails of solitons at n — > +00, we 
obtain the corresponding conditions, 

Cij-\C<2,j = CzjCtj-i — G'2jC2j-i — G'2j-iC2j — O, j = l,2,..., N, 
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on the residue matrices. For instance, the one-soliton solution of ( |13|) is computed as 
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In the continuum limit of time (St — > 0) , the one-soliton solution reduces to that of the semi-discrete 
CNLS equations |J due to the relations, 

a = St sinh 2W cos 29 + 0(8t 2 ), 
b = St cosh2W sm26 + 0(St 2 ). 

It is noteworthy that we can obtain the most general soliton solutions of the model since we have 
employed the ISM for the analysis. 

Conserved quantities of the full-discrete CNLS equations ( p"2| ) under the decaying boundary 
conditions can be computed recursively by expanding the time-independent quantity, logdet A(z), 
with respect to 1/z. Here we list first three conserved densities, 

iog(i-E?i J M J) ), 

3 

cos2aE(9i+i^' ) +^'V« 1 ) + isin2 a E('?n+i^ ) " <^'M*+i), 



3 j 3 

- i{cos2aE(^i^' ) + 1^1) + isin2aE(^i^' ) " <^i)} 2 



/ j In+l'n+l / y In ' n ■ 

3 j 



8 



4 Concluding remarks 



We have investigated an integrable full-discretization of the coupled nonlinear Schrodinger (CNLS) 
equations from the point of view of the inverse scattering method (ISM) for the first time in this 
article. We have found a Lax pair for full-discrete CNLS equations and performed the ISM under 
appropriate conditions with the help of the transformation (|l6|). As a result, we can solve the 
initial- value problem, and obtain soliton solutions and conserved quantities. 

There are plural schemes for integrable full-discretization of one-component nonlinear Schrodinger 
(NLS) equation. We have some freedom to determine the linearized dispersion relation of full- 
discrete NLS equation (see ||). However, as far as we have considered, there is no such freedom 
in the multi-component case, which leads us to ( |T^ ) as full-discrete CNLS equations. It is to be 
noted that ([l2|) can be cast into a few types of full-discrete CNLS equations by coordinate trans- 
formations and redefinitions of the auxiliary variable j n (see W for the work of Hirota, Ohta and 
Tsujimoto on full-discrete NLS equation). 

Finally, we comment that discrete CNLS equations have been studied by using the direct 
method p, 0, |l0) ■ In (|] , Hirota proposed the full-discrete coupled modified KdV equations ( |l8| ) 
and gave a (1 + l)-soliton solution in the two-component case. Ohta obtained an ./V-soliton solution 
for the semi-discrete CNLS equations in terms of Pfaffian ||. However, the ./V-soliton solution 
is not a general one since the initial soliton polarizations in the solution are either parallel or 
orthogonal. More recently, by using another Pfaffian representation, Ohta proposed more general 
soliton solutions without a restriction on soliton polarizations in the full-discrete case . He also 
obtained a Lax pair of different form in |Io| by means of an alternative approach. 
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